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LONG-TERM  GOALS 

Providing  particle  characteristics  such  as  particle  size  distribution  (PSD)  and  particle  composition  (n^) 
from  ocean  color  data  greatly  facilitates  classifying  water  types  and  predicting  optical  properties  of 
natural  waters.  This  information  impacts  many  Navy  warfare  areas  which  depend  on  active  and 
passive  imagery  or  that  require  knowledge  of  particle  types  and  their  dispersion.  Potential  methods  to 
predict  particle  characteristics  using  the  inherent  optical  properties  (lOPs)  has  been  proposed  (Boss  et 
al.  2001;  Twardowski  et  al.  2001).  Haltrin  and  Arnone  (2003)  put  forward  an  algorithm  to  separate 
particle  type  into  large  and  small  particles  and  terrigenic  versus  organic  source  and  concentration  of 
particles  using  remote  sensing  data  using  spectral  particulate  coefficient  information.  This  work 
however  relied  heavily  on  assumptions  between  the  backscattering  to  scattering  ratio  and  the  ability  to 
iterate  using  scattering  efficiencies  to  derive  total  particle  number  given  a  Junge  slope  via  Twardowski 
et  al.  The  ultimate  goal  of  the  work  reported  here  is  to  examine  these  algorithms  and  their  assumptions 
in  order  to  provide  more  confidence  in  the  methods  used  to  retrieve  in-situ  particle  and  optical 
properties  from  either  passive  or  active  remote  sensing  and  from  either  underwater,  airborne,  or 
satellite  platforms.  Knowing  the  optical  and  particle  properties  of  the  water  can  improve  the  predicted 
performance  of  active  and  passive  imaging  systems  and  our  ability  to  differentiate  possible  biological 
threats  in  an  area.  It  is  the  goal  of  this  research  to  strive  to  produce  such  data  that  will  aid  in  the 
characterization  of  the  battlespace  environment  for  electro-optical  systems  deployed  by  the  Navy  for 
Mine  Counter  Measures,  Naval  Special  Warfare,  Intelligence  Surveillance  and  Reconnaissance,  and 
Anti-Submarine  Warfare. 

OBJECTIVES 

To  improve  our  estimates  of  coastal  optical  properties  from  remote  sensing  such  a  capability  should 
take  advantage  of  the  influence  of  the  PSD  and  bulk  Up  on  the  inherent  optical  properties  (lOPs)  of  the 
water  column.  The  objective  of  this  effort  is  to  determine  the  variability  in  the  relationship  between 
measured  lOPs  and  various  measured  particle  characteristics  for  lake  and  coastal  waters.  In  order  to 
use  relationships  such  as  those  of  Boss  et  al.  (2001)  and  Twardowski  et  al.  (2001)  there  must  be 
confidence  in  the  Junge  functions  description  of  coastal  environments  and  the  impact  of  composition 
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on  scattering  and  attenuation.  The  objeetives  for  this  effort  are:  1)  to  determine  the  robustness  of 
existing  models  to  prediet  optieal  and  partiele  eharaeteristies,  2)  determine  whether  models  that 
estimate  refraetive  index  and  shape  properties  ean  be  used  with  partiele  eoneentration  to  derive 
observed  in-situ  optieal  properties  (scattering  and  baekseattering),  and  3)  identify  where  present 
models  should  be  adjusted  to  better  retrieve  particle  eharaeteristies.  The  work  will  include  attempts  of 
elosure  between  partiele  size  distribution,  partiele  eomposition,  ealeulated  baekseattering  and  total 
scattering,  and  measurements  of  the  volume  seattering  funetion.  The  measured  volume  seattering 
funetion  for  various  water  types  are  necessary  to  evaluate  disparities  in  using  remote  sensing 
refleetanees  to  retrieve  optieal  properties. 

APPROACH 

The  New  York  Finger  Lakes  are  used  as  a  test  area  in  order  to  cover  a  wide  dynamic  range  of  inherent 
optieal  properties,  partiele  eharaeteristies,  and  possible  remote  sensing  differenees.  These  lakes  range 
from  meso-oligotrophie  to  hyper-eutrophie  and  from  ehlorophyll  dominated  to  inorganie  mineral 
dominated  particle  types.  The  main  approaeh  is  to  measure  the  partiele  characteristics  and  the  in-situ 
optieal  properties  for  these  various  water  types.  Lakes  were  ehosen  sinee  being  enelosed  they  ean  be 
largely  eonstrained  and  represent  numerous  partiele  types  due  to  fluvial  input,  plankton  diversity,  and 
varying  eoneentrations  of  organie  and  minerogenie  partieles.  Even  with  these  diverse  eharaeteristies 
they  are  for  the  most  part  stable  over  sampling  intervals.  Lakes  are  largely  impacted  by  “episodic” 
events  and  unlike  eoastal  oeeans  the  physieal  foreing  is  mostly  point  souree  inputs,  stratifieation,  and 
surfaee  mixed  layers. 

Lake  sampling  was  conducted  from  27  July  to  4  August  2005  during  a  period  when  the  Finger  Lakes 
are  dominated  by  different  partiele  types  (inorganie  CaCOs,  phytoplankton  (diatoms  versus  bluegreen), 
and  elay).  The  lakes  inelude  Onondaga  Lake,  Otiseo  Lake,  Owaseo  Lake,  Skaneateles  Lake,  Seneea 
River,  and  Cross  Lake.  Seneea  River  stations  represent  areas  where  phytoplankton  ehanges  and 
partiele  eharaeteristies  (inorganie  and  organie)  ehange  rapidly  due  to  zebra  mussel  presenee  with  high 
seattering  b  (b(650)>7  m'^).  In  eontrast  Skaneateles  is  a  low  produetion  lake  dominated  by 
phytoplankton  with  little  inorganie  and  low  seattering  (b(650)<0.5  m'^).  The  optieal  measurements 
include  1)  laboratory  spectral  absorption,  (absorption  by  detritus,  CDOM,  and  phytoplankton),  2)  in- 
situ  speetral  total  absorption,  fdtered  absorption,  and  total  attenuation  (Wetlabs  Ine.  aes),  3)  in-situ 
speetral  scattering  measurements  at  seleeted  angles  (WetLabs  Ine  ECO-VSF,  and  9-ehannel 
baekseattering  sensor),  4)  in-water  speetral  downwelling  irradianee  and  upwelling  radianee  (Satlantie 
Ine.,  Hyperpro),  5)  laboratory  measurements  of  the  volume  scattering  function  at  8  wavelengths 
(MVSM,  NRLSSC),  and  6)  laboratory  and  in-situ  partiele  size  distribution,  and  laboratory 
measurements  of  ehemieal  eomposition.  Exeept  for  Onondaga  Fake,  samples  were  taken  at  the 
surfaee,  for  Onondaga  Fake  samples  were  eollected  at  the  surfaee  and  also  near  the  bottom  (8  m). 

Particle  eharacterization  was  earried  out  on  bulk  samples  in-situ  and  on  laboratory  samples.  Partiele 
properties  measured  inelude  the  partiele  size  distribution  (Sequoia  Seientifie  FISST),  total  and  volatile 
suspended  sediment  eoneentrations,  and  ehemieal  eonstituents  of  the  partieles.  Properties  of  the 
partieles  inelude  measurements  of  the  partiele  size  distribution,  inorganie  and  organie  suspended 
material,  and  the  ehemieal  composition  and  a  shape  index.  The  ehemieal  composition  of  surface 
samples  that  were  measured  with  the  MVSM  were  also  analyzed  using  a  Seanning  Eleetron 
Mieroseope  modified  with  an  X-ray  energy  speetrometer  (SAX).  This  instrument  provides  partiele 
size,  shape  index,  and  ehemieal  composition  (but  has  been  demonstrating  some  bias  with  respeet  to 
inorganie  and  organie  particles).  The  ehemieal  eomposition  information  is  being  used  with  partiele 


size  data  to  estimate  the  bulk  relative  refractive  index  of  the  sample  as  well  as  the  scattering  and 
backscattering  via  Mie  theory.  The  particle  size  distribution  (PSD)  was  measured  with  the  SAX  and 
then  the  slope  derived  over  the  1 .25  um  and  higher  range  fitting  to  the  Junge  functional  form.  This 
was  necessary  due  to  the  “non-Junge-like”  distribution  as  will  be  presented  below.  The  spectral  slope 
of  the  attenuation  coefficient  as  discussed  in  Boss  et  al.  (2001)  and  Twardowski  et  al.  (2001)  is 
compared  with  the  particle  size  distribution  slope.  The  results  are  then  used  to  help  interpret  data 
collected  using  the  MVSM  in  the  backward  direction. 

WORK  COMPLETED 

Water  samples  have  been  collected  from  ten  locations  in  New  York;  1)  Owasco  Lake,  2)  Skaneateles 
Lake,  3)  Otisco  Lake  (North  basin)  4)  Otisco  Lake  (South  basin),  5)  Onondaga  Lake  (1  and  8  meter 
depths),  6)  Seneca  River  Buoy  424,  7)  Cross  Lake  (North  area),  8)  Cross  Lake  (South  area),  9)  Seneca 
River  Buoy  409,  and  10)  Seneca  River  Buoy  B317.  For  each  station  laboratory  measurements  of 
particulate  and  dissolved  absorption  using  a  spectrometer,  as  well  as,  particle  size,  asphericity  index, 
and  chemical  composition  using  the  SAX  have  been  made.  MVSM  scattering  function  was  measured 
at  nearly  all  the  sites.  In-situ  data  collected  includes  the  total  attenuation,  particulate  absorption, 
dissolved  absorption,  particle  size  distribution  (via  LISST),  and  spectral  backscattering.  In  these  turbid 
waters  (c  ranging  from  <1  to  >  8  m'^)  there  was  uncertainty  in  values  of  “b”.  To  help  resolve  this  a 
comparison  between  laboratory  (fdter  pad,  and  dissolved)  and  field  measured  (Wetlabs  acs  with  and 
without  filter)  absorption  was  conducted.  While  some  variance  was  present  in  the  blue  end  of  the 
spectra,  overall  this  gave  confidence  to  the  SAX  and  spectral  Cp  comparison.  In  addition 
measurements  of  the  volume  scattering  function  (MVSM)  were  made  since  multiple  scattering  issues 
have  been  previously  suggested  in  this  instrument  but  such  effects  are  likely  in  the  other  in-situ 
instrumentation  as  well. 

RESULTS 

A  large  range  in  optical  conditions  was  observed;  which  was  a  primary  reason  for  the  selection  of  the 
Finger  Lakes  and  Seneca  River  for  some  of  the  inherent  optical  property  and  particle  characteristics 
investigation.  Figure  1  shows  that  particulate  attenuation  at  550  nm  ranges  from  0.52  in  Skaneateles  to 
9.0  in  Seneca  River.  Within  the  Seneca  River  stations  there  was  an  eightfold  difference  between  the 
stations  impacted  by  the  zebra  mussels.  There  is  also  clear  evidence  in  figure  1  of  the  difference  in 
particle  types  as  demonstrated  by  the  slope  in  the  particulate  attenuation  and  the  SAX  chemical 
analysis.  SAX  data  showed  clay  to  be  the  largest  single  component  of  the  minerogenic  fraction  in  all 
of  the  sites  reaching  over  70%  at  the  buoy  station  424  on  Seneca  River.  Calcium  Carbonate  was  a 
dominate  contributor  in  Otisco  Lake,  while  diatoms  dominated  much  of  the  projected  area  per  unit 
volume  of  high  refractive  index  material  in  Skaneateles. 

Twardowski  (2001)  and  others  have  suggested  that  the  slope  of  the  particulate  attenuation  can  be  used 
to  describe  the  particle  size  distribution.  As  figure  1  demonstrates  a  large  variability  in  this  slope  is 
present  in  these  test  areas.  Figure  2  represents  an  analysis  of  particle  size  distributions  for  these 
various  water  bodies  over  the  range  of  0.2  to  30  pm;  the  best  operating  range  for  the  SAX 
instrumentation.  Some  loss  of  organic  contributions  is  expected  with  the  SAX  but  the  overall  particle 
size  distribution  (PSD)  is  not  expected  to  change  significantly  in  these  waters.  Note  that  these  waters 
do  not  follow  Junge  type  distribution  but  rather  the  2  component  shape  of  Risovic  (1993).  In  order  to 
best  compare  the  spectral  cp  slope  and  the  PSD  slopes,  a  Junge  type  distribution  was  fit  over  the  region 
of  2  to  30  pm.  This  still  left  the  slope  of  the  PSD  below  a  value  of  3  (range  of  2.5  to  2.9  for  these 


stations)  and  outside  the  normal  oeeanographie  range.  This  is  not  surprising  given  the  non-Junge  like 
nature  of  the  PSD  from  the  SAX. 
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Figure  1.  Particulate  attenuation  for  each  study  site  [graph:  ac-s  derived  particle 
attenuation  data  showing  changing  slope  in  particulate  attenuation,  panel  A 
representing  the  low  attenuation  waters  (note  very  low  slope  for  Skaneateles  Lake)  and 

panel  B  the  higher  attenuating  waters] 
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Particle  Size  Distribution  from  SAX 
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Figure  2.  Particle  size  distribution  from  SAX  [graph:  the  particle  size  distribution  over 
range  of  0.2  to  30  fim  for  each  sample  showing  a  non-Junge  functionality  with  a  peak 

in  the  particle  size  at  about  4  to  5  /im] 


To  determine  the  ability  of  the  SAX  to  retrieve  the  baekscattering  and  scattering  coefficient  the  particle 
size  distribution  and  the  separation  of  inorganic  and  minerogenic  particles  was  used.  From  the 
chemical  composition  a  bulk  relative  refractive  index  was  calculated.  These  values  compared  well 
with  those  determined  using  particulate  baekscattering  to  total  scattering  ratios  from  the  in-situ  data. 
The  range  of  values  was  from  1.07  for  Skaneateles  (phytoplankton  dominated)  to  1.16  in  Otisco  North 
where  calcium  carbonate  persisted.  Most  other  stations  were  about  1.14  to  1.15.  In  contrast  using  the 
baekscattering  to  total  scattering  ratio  the  minimum  refractive  index  was  1.10  for  Skaneateles 
increasing  to  1.19  in  the  Seneca  River. 


Figure  3.  Comparison  of  scattering  and  baekscattering  at  650  nm  measured  in  situ 
versus  that  calculated  from  particle  characteristics  for  minerogenic  fraction  for  a 
series  of  New  York  Finger  Lakes  [graphs:  measured  particle  scattering  b  and 
baekscattering  bbp  versus  that  calculated  from  SAX  based  on  Mie  Theory  using  SAX 
size  and  chemical  composition,  panel  (A)  measured  b  and  SAX  bm  with  inset  showing 
adjustments  to  scattering  based  on  asphericity  factor  (ASP),  (B)  measured  bbp,  and  bbm 
by  SAX  with  average  asphericity  factor,  and  panel  (C)  b  bp/bp  measured,  and  bbm/bm  by 
SAX  (m=  mineral  component,  o=organic),J 


However  the  SAX  particle  size  distribution  and  composition  data  was  also  used  to  calculate  the  total 
scattering  and  backscattering  coefficients  at  650  nm  (Figure  3)  via  Mie  Theory.  When  this  was  done 
the  results  showed  that  while  the  SAX  inversion  retrieved  much  of  the  particulate  backscattering,  there 
was  a  large  discrepancy  in  the  total  scattering  (and  therefore  the  ratio  of  backscattering  to  total 
scattering).  An  asphericity  factor  (ASP)  was  calculated  from  the  SAX  data  and  represents  a  measure 
of  deviation  from  a  sphere.  This  did  not  impact  the  results  of  either  the  scattering  or  backscattering 
significantly.  For  the  most  part  the  difference  is  likely  due  to  the  inability  of  the  SAX  to  properly 
account  for  all  of  the  organic  fractions  whose  size  classes  influenced  forward  scattering  significantly 
(see  panel  A).  This  is  expected  and  also  demonstrated  in  the  “SAX-minerogenic”  fraction  which 
shows  a  particulate  hb/h  ratio  more  closely  aligned  with  values  expected  for  highly  refractive 
compounds. 

The  MVSM  measurements  often  showed  a  distinct  feature  in  scattering  in  the  backward  direction  from 
140  to  170  degrees  (Figure  4).  In  order  to  evaluate  this  change  in  shape  of  the  scattering  function  the 
measurements  were  normalized  to  total  scattering  b,  and  to  scattering  at  120  degrees.  This  was  to 
remove  spectral  dependencies  that  were  evident  in  the  data  as  well  as  give  a  “hinge”  point  that  is  close 
to  that  often  used  in  other  instrumentation  for  determining  total  backscattering. 
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Figure  4.  Onondaga  surface  and  8m  backscattering  function 
[graph:  between  130  and  1 70  degrees  there  is  a  marked  difference  in  the  normalized 
scattering  function  between  the  surface  and  subsurface  samples] 

This  change  in  the  shape  of  the  Volume  Scattering  Function  was  thought  in  previous  years  to  be  fully  a 
function  of  instrument  artifacts.  However  within  this  past  year  with  careful  studies  there  is  now  some 
evidence  to  suggest  that  this  difference  may  in  fact  represent  a  convolution  of  instrument  sensitivities 
and  actual  differences  in  the  shape  function.  This  feature  is  being  addressed  more  fully. 


IMPACT/APPLICATIONS 


The  impact  of  the  work  is  that  closure  between  particle  optics  and  bulk  optical  properties  will  provide 
a  direct  tie  to  use  of  remote  sensing  for  characterizing  coastal  environments.  The  research  also 
demonstrates  how  robust  the  methods  used  to  derive  backscattering  for  various  particulate  assemblages 
are  in  lakes.  While  seemingly  not  applicable  to  marine  dynamics,  many  coastal  areas  or  semi-enclosed 
basins  resemble  in  their  particle  types  and  dynamics.  The  use  of  chemical  composition  and  particle 
properties  to  evaluate  bulk  properties  will  impact  how  we  interpret  remote  sensing  signals,  particularly 
in  areas  of  denied  access  or  for  areas  where  potential  contamination  is  suggested.  Lastly  for  Navy 
applications  the  volume  scattering  variability  and  the  impact  of  particle  types  are  critical  for 
determining  blurring  properties  that  will  be  used  for  system  performance  modeling  of  Navy  imaging 
systems.  If  particle  types  can  be  deduced  from  remote  sensing,  then  this  information  can  be  utilized 
for  defining  biological  growth  areas,  sediment  loading,  and  points  of  discharge  in  satellite  imagery. 

TRANSITIONS 

The  volume  scattering  data  has  been  sent  to  the  Canadian  Defense  Research  Development  Center  (Dr. 
George  Fournier)  for  application  to  the  LUCIE  system  for  analysis  on  the  impact  to  blurring  on  this 
system  of  particles  with  different  characteristics;  in  addition,  the  non-Junge  scattering  functions  impact 
on  blurring  is  being  evaluated. 

The  combined  particle  analysis  and  lake  inherent  optical  properties  is  being  used  by  the  Upstate 
Freshwater  Institute  (Dr.  S.  Effler)  for  long-term  monitoring  of  New  York  Finger  Fakes. 

RELATED  PROJECTS 

Improved  sensor  performance  using  ‘through  the  sensor’  techniques  (ONR/NRL  6.2  POC-A. 
Weidemann)-  goal  is  to  derive  as  much  information  about  the  volume  scattering  function  and/or 
particle  types  from  gated  laser  system  return  at  multiple  distances  in  order  to  improve  image  quality  in 
subsequent  images 
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